Universal scaling laws of chaotic escape in dissipative multistable systems subjected to 

autoresonant excitations 
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A theory concerning the emergence and control of chaotic escape from a potential well by means of 
autoresonant excitations is presented in the context of generic, dissipative, and multistable systems. 
Universal scaling laws relating both the onset and lifetime of transient chaos with the parameters of 
autoresonant excitations are derived theoretically using vibrational mechanics, Melnikov analysis, 
and energy-based autoresonance theory. Numerical experiments show that these scaling laws are 
robust against both the presence of noise and driving re-shaping. 
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Introduction.— ¥iiic&^e from a potential well is an old 
problem with wide-ranging implications where the in- 
terplay of noise, dissipation, deterministic driving, and 
quantum uncertainty has given rise to unexpected and 
intriguing phenomena such as coherent destruction of 
tunnehng [1] and stochastic resonance [2]. Other spe- 
cific examples are pulse-shape-controlled tunneling [3], 
shot-noise-driven escape in Josephson junctions [4], and 
thermally induced escape [5]. While most previous in- 
vestigations on driven escape have restricted themselves 
to purely periodic drivings, chirped excitations also have 
a demonstrated effectiveness. Chirped lasers, for exam- 
ple, can reduce the intensity required for infrared multi- 
photon dissociation of diatomic molecules to an experi- 
mentally realizable intensity range [6]. And chirped op- 
tical pulses have been shown to enhance charge flow in 
molecular-tunneling junctions [7]. In spite of the impor- 
tance of chirped excitations [8] , to the best of the author's 
knowledge, laws governing (some of) the associated es- 
cape scenarios remain to be revealed. The key to un- 
derstanding the aforementioned effectiveness of chirped 
excitations is essentially the autoresonance (AR) mech- 
anism: AR induced energy amplification in nonlinear, 
driven, and deterministic systems occurs when the sys- 
tem continuously adjusts its amplitude so that its in- 
stantaneous nonlinear period matches the instantaneous 
driving period of the chirped excitation. Initially studied 
in the context of a Hamiltonian description, AR phe- 
nomena have been well known for about half a century 
and have been observed in particle accelerators, plane- 
tary dynamics, atomic and molecular physics, and non- 
linear oscillators [9], to cite a few examples. Regarding 
dissipative systems, an energy-based AR (EBAR) the- 
ory has recently been proposed and applied to the case 
where the system crosses a separatrix associated with its 
underlying integrable counterpart [10]. Since in such an 
escape situation the appearance of transient chaos asso- 
ciated with the occurrence of homoclinic bifurcations is 
an ubiquitous phenomenon, the question naturally arises: 
How does AR control the chaotic escape scenario, i.e., the 
onset and lifetime of transient chaos in generic dissipative 
multistable systems? 



Theory.— In this Letter, this fundamental problem is 
studied in the context of the family of systems 

x + dU/dx = -5x + 7 sin [Vl {t) t] , (1) 

where U{x) is a generic multistable potential (see Fig. 1), 
and Q.{t) = u: + ant"", n = 1,2,..., is a time-dependent 
frequency with being the n-th-order sweep rate. One 
expects [10] 7 sin [fi (t) t] to behave as an effective au- 
toresonant excitation whenever the sweep rate is suffi- 
ciently low (adiabatic regime), which is the case con- 
sidered throughout the present paper in order to apply 
Melnikov analysis (MA) to autoresonant excitations [11]. 
Also, the damping and autoresonant excitation terms are 
taken to be small amplitude perturbations of the under- 
lying integrable system so as to deduce analytical expres- 
sions for the scaling laws relating both the onset time ti 
and lifetime t of transient chaos with the parameters of 
the AR excitation from the MA [12] results. The appli- 
cation of MA to any homoclinic (or heteroclinic) orbit 
[xh,j (t) , Xhj (t)] of the unperturbed (5 = 7 = 0) coun- 
terpart of Eq. (1) involves calculating the corresponding 
Melnikov function (MF) 

(to) =-D + jRo (iJ, to) + la^Ri {uj, to) + O {jal) , 

(2) 

where D = S J^^Xf^ ,^ (t) dt > 0, Ro{uj,to) = 
J^^Xhj {t)sm[uj{t + to)]dt, Ri{uj,to) = 

+ ^0)" (0 cos [w (t -I- to)] f^^- For a purely 
periodic excitation (a„ =0), it is assumed in the follow- 
ing that the system does not exhibit transient chaos for 
a given set of parameter values, i.e., the MF (2) has no 
simple zeros: 



-D > 7 



max Ro (o^, tn] 

to6 K 



7^0 (t^) , 



(3) 



where Ao (lo) represents a chaotic-threshold function 
whose generic behaviour is shown in Fig. 1. Indeed, 
homoclinic chaos is not possible for sufficiently small 
Lu since a purely harmonic excitation becomes a con- 
stant when cj — > 0. For sufficiently high uj, the dy- 
namics can be analyzed using the vibrational mechan- 
ics approach [13] by separating x{t) = z{t) + ipit), 
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FIG. 1: (color online) Generic multistable potential U{x) vs 
X, and energy levels corresponding to different separatrices 
(dashed lines). The inset shows a generic chaotic threshold 
function Ao (uj) vs u> and the range [u}th,ij'^th,2] in which ho- 
moclinic chaos is expected [see Eq. (3)]. 



where z{t) represents the slow dynamics while tp {t) is 
the fast oscillating term: ip (t) = V'o cos {u)t + <Pq) with 

■(/,'() = \/ (5^ + Lo"^ / i^ujS"^ — uj^) , = arctan ((5/a^')- 
averaging out ip (t) over time, the slow reduced dynamics 
of the system becomes 

z + dV/dz = —Sz, 

dV/dz = T-^ I g[z + 'il}QCOs{ujt + ipQ)]dt,(A) 
Jo 

where g (x) = dU{x)/dx,T = 2tt/uj, i.e., that of a 
purely damped system, and hence homoclinic chaos is 
not possible when uj oo. Note that this demon- 
strates that equilibria are the only attractors of the sys- 
tem for «„ > 0. Thus, one concludes that the proper- 
ties Ai) {iu 0, oc') = 0, Aq (lu) ^ imply via the ex- 
treme value theorem (Weierstrass' theorem [14]) that 
the generic chaotic-threshold function (w) presents at 
least one maximum (the case shown in the inset of Fig. 1). 
Equation (2) indicates that the MF Mj (to) has simple ze- 
ros (at sufficiently large values of to) for any positive value 
of the sweep rate because of the factor (t -I- ^o)" appearing 
in the definition of Ri {u!,to). Physically, this means that 
after a sufficiently long time, ti, which depends upon the 
system's parameters and initial conditions, the instanta- 
neous frequency of the autoresonant excitation reaches 
a threshold value Q (t,;) = LUih.i that is the lowest fre- 
quency satisfying the relationship D/^ = Aq {uJth,i) (see 
Fig. 1, inset), i.e., the threshold condition for the onset of 
chaotic behaviour. It is worth noting that this occurs for 
any initial conditions because of the AR- induced increase 
of the system's energy. Thus, the condition Cl (ti) = Uth,i 
implies the scaling law 



1/n 



(5) 



for the onset time of transient chaos. Similarly, the life- 
time r of the chaotic transients can be estimated from 



the instantaneous frequency [ti + t) ^ that is 

the lowest frequency satisfying the relationships D/^ = 
{u)th,2) ,iOth,2 > wth.i (see Fig. 1, inset), i.e.. 



T ~ 



(Wt/i,2 - W)^^"" - {(jJth,! - W) 



1/n 



a. 



-1/n 



(6) 



We see that the universal scaling laws (5) and (6) 
are inverse-power laws containing the two parameters 
{w,an) that control the autoresonant excitation while 
the critical exponent is the inverse of the chirp order n. 
Notably, these scaling laws also contain the dependence 
upon the particular potential, initial potential well, and 
dissipation and excitation strengths through the thresh- 
old frequencies ujth,i,i^th,2- Since the mid 1980's, critical 
exponents of chaotic transients have been discussed the- 
oretically in the context of crises in dissipative maps [15]. 
The present theory establishes that an inverse-power law 
remains valid also for dissipative multistable flows sub- 
jected to autoresonant excitations. The connection be- 
tween such results for maps and the present ones can be 
readily understood by assuming that for a purely periodic 
excitation (q!„ = 0) there exists a chaotic attractor for a 
certain frequency {u)th,i < OJ < L0th,2) instead of a peri- 
odic attractor {oJth,i > w) (see Fig. 1, inset). According 
to the above discussion, the chaotic attractor disappears 
for an > via a boundary crisis and chaotic transients 
appear instead, with a = Q.n,c = being the critical value 
for the crisis. Therefore, since the non-existence of sim- 
ple zeros of the MF is a sufficient condition for the disap- 
pearance of transient chaos, the lifetime of these chaotic 
transients follows the universal inverse-power law 



a/»i -i/u 



(7) 



Now, applying EBAR theory one can deduce scaling 
laws relating both the onset time and lifetime with the 
autoresonant excitation amplitude. In particular, for 
Duffing-like potentials one has that the optimal ampli- 
tude scales as 7 ~ [3" (n + 1)!]'/^'"+'^ [10] 
and hence the scaling law (7) for example becomes 

T ~ [(a;tfc,2-w)(n+l)!]'/"7-(2"+2)/(3"\ where one 
sees that the critical exponent ranges from 4/3 for a lin- 
ear chirp to 2/3 for the highest-order chirps. 

Robustness vs noise and re-s/iapm^.— Numerical 
experiments confirmed the accuracy of the above scal- 
ing laws in different systems. Representative results 
corresponding to a dimensionless Duffing oscillator 
X ^ X - x^ - Sx + jsn [2K{ni)n (t) t/n; m] + aN (0, 1) 
are shown in Figs. 2-5 for illustrative purposes. Here, 
sn(-;m) is the Jacobian elliptic function of parameter 
m E [0, 1] with K{m) being the complete elliptic integral 
of the first kind. One has sn(-;m = 0) = sin(-) while, 
in the other limit, sn(-;m = l) is the square- wave 
function [16]. Also, a N is a white-noise term 

with N (0, 1) being a random variable with Gaussian 
distribution, zero mean, and variance 1, while cr > 
controls the noise strength. For the purely deterministic 
(ct = 0) case of a sinusoidal (m = 0) excitation with a 
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linear chirp {n = 1), for example, and after calculating 
the resulting integrals by residues, one obtains the 
corresponding MF: M^^^.^.^^ (to) = - D ^ Bq cos (uto) T 
ai [Bito cos (Luta) - {B2 + B^tl) sin (oito)] + O (70;?) , 
where the sign + (— ) corresponds to the 
right (left) homoclinic orbit of the unper- 
turbed Duffing oscillator ((5 = 7 = 0), and 
D = 4(5/3, Bq = V27r7Wsech(7ru;/2), Bi = 
4\/27 [{tt/2) - {■k'^uj/^) tanh(7ra;/2)] sech(7ra;/2), B2 = 
\/27 (ti'^/S) scch^ (7r(jj/2) [4sinh {ttlo) — ttuj cosh (ttuj) + Sttlo]. 
The results for the case of a periodic (chaotic) attractor 



FIG. 3: (color online) Autoresonant response of a dimension- 
less Duffing oscillator for a linear chirp ai = 3 x 10"'* (see the 
text), (a) Position vs time and final equilibrium (dashed line), 
(b) Energy and average energy (over a few periods 2-k/lj, thick 
black line) vs time, (c) Phase space trajectory and period-1 
attractor existing at a„ = (thick black line), (d) Energy 
vs position (dashed line indicates the separatrfx energy level 
as in version (b)). Remaining parameters are the same as in 
Fig. 2. 



FIG. 2: (color online) Onset time (circles), ti, and life- 
time (squares, inset), r, of transient chaos corresponding to 
a dimcnsionless Duffing oscillator (sec the text) for m = 
0, a = 0, 5 = 0.5, 7 = 0.4, and oj = 0.493. Hence, ujth,2 = 
1.11736, wtfe,i = 0.49371 from M±^^^.„g (to) (see the text). 
Solid lines indicate fits according to the scaling laws (5) and 
(6), respectively. 

existing at q;„ = are shown in Fig. 2 (Fig. 4). While a 
linear chirp is considered in Figs. 2 and 3, Fig. 4 shows 
the results for quadratic and cubic chirps. Figures 2 and 
4 show excellent agreement between numerical results 
and the predicted scaling laws for both the onset time 
and the lifetime of transient chaos. Figure 3(b) shows 
the AR-induced increase of energy over time until the 
separatrix energy level is reached, indicating the onset 
of transient chaos (see Figs. 3(a) and 3(d)), followed 
by a decrease of energy when the dynamics becomes 
effectively purely dissipative (see Fig. 3(c)). Note that 
the average energy remains roughly constant during 
transient chaos. The robustness of the scaling laws 
versus both the presence of additive noise and rc-shaping 
of the autoresonant excitation is shown in Fig. 5. Note 
that, for the elliptic excitation, u!th.i,i — 1,2, are 
functions of m, which is the reason for the different 
prefactors in the fits for m = 0.995 and m = 1 — 10~^^ 
(Fig. 5, bottom). 

Conclusions.— In sum, universal inverse-power laws re- 
lating both the onset and lifetime of transient chaos and 



FIG. 4: (color online) Lifetime of transient chaos correspond- 
ing to a dimcnsionless Duffing oscillator (see the text) for 
quadratic (n = 2, squares) and cubic (n = 3, circles) chirps. 
Solid lines indicate fits according to the scaling law (7). Also 
shown is the chaotic attractor existing at 02,3 = 0. Fixed pa- 
rameters: m = 0, (T = 0, (5 = 0.154,7 = 0.2, u; = 1.1. Hence, 
0Jth,2 = 1.71404 from M^^ , {to) (see the text). 
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the parameters of escape-inducing autoresonant excita- 
tions have been theoretically derived for generic, dis- 
sipative, and multistable flows. Numerical simulations 
showed the robustness of these scaling laws against both 
the presence of additive noise and driving re-shaping. 
Since the critical exponents were foimd to solely depend 
on the chirp's order, such scaling laws are expected to re- 
main valid for even more general dissipativc systems in- 
cluding spatiotemporal chaotic systems [17]. The present 
results can be readily tested experimentally (for example 
in mechanical and laser systems), and can find applica- 
tion to optimally control elementary dynamic processes 
characterized by chaotic escapes from a potential well, 
such as diverse atomic and molecular processes, trans- 
port phenomena in dissipative lattices, and control of the 
electron dynamics in quantum solid-state devices. 



FIG. 5: (color online) Lifetime of transient chaos correspond- 
ing to a dimensionless Duffing oscillator (see the text) in the 
presence of noise (top panel, m = 0, 7 = 0.4) and subjected 
to elliptic excitations (bottom panel, cr = 0, 7 = 0.3). Solid 
lines indicate fits according to the scaling law (6). Fixed pa- 
rameters: 6 = 0.5, ui = 0.493. 
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